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Evolution of the X-ray luminosity and metallicity of 
starburst blown superbubbles 
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ABSTRACT 

We calculate the time-dependent metal production expected from starbursts and use 
them as boundary conditions in our 2D simulations of evolving superbubbles. We as- 
sume that the produced metals (oxygen and iron) thoroughly mix with the ejected 
stellar envelopes, and/or with the matter thermally evaporated from the superbubble 
cold outer shell. The metal production process determines the time-dependent metal- 
licity in hot superbubble interiors, and leads to values of Z > Z Q when oxygen is 
used as tracer, and under-solar when the metallicity is measured with respect to iron. 
In either case, the enhanced metallicity boosts the X-ray emissivity of superbubbles, 
bringing theory and observations closer together. 

Key words: ISM: abundances, ISM: bubbles, ISM: hydrodynamics, galaxies: star- 
bursts, X-rays: bursts 
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1 INTRODUCTION 

The formation and evolution of galaxies is one of the most 
captivating problems of modern astrophysics. With present- 
day techniques, young star clusters can be observed at large 
look-back times and much work has been devoted to the de- 
tailed studies of the optical and near UV properties of such 
systems at intermediate and high redshifts. The defining 
characteristic of starburst galaxies is their spectrum which 
is the emission lines from HII regions in extremely young 
(< 10 Myr) star forming regions. While the UV radiation 
from the massive stars excite and ionize the associated HII 
regions, the stellar winds and supernova explosions lead to 
giant hot bubbles in the interstellar medium (ISM). 

Most of the information about the properties of the 
"warm ISM" in young giant star-forming regions comes from 
the analysis of their emission lines. Such observations in- 
dicate that the warm ISM in starburst galaxies is not yet 
contaminated by the metals ejected by the present burst. 
However, little is known about the "hot ISM", where most 
of the newly produced metals are presumably located. Data 
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from ROSAT, ASCA and BeppoSAX , although inconclu- 
sive, seem to suggest sub-solar iron abundance and about so- 
lar abundances for the a-elements (Bauer & Bregman 1996, 
Ptak et al. 1999, Persic et al. 1998, but see Dahlem, Weaver 
& Heckman 1998 and Weaver et al. 1999). However, The 
abundances determined from the optical emission lines, in 
some of the best studied galaxies, is considerably higher than 
those derived from the X-ray data. This conflict is likely to 
be resolved soon, with the coming new high resolution spec- 
troscopic observations of Chandra and XMM. 

The mechanical energy from evolving starbursts is 
known to lead to the formation of superbubbles in the ISM. 
This happens as the violently ejected matter from winds and 
supernovae becomes thermalized at a reverse shock. This 
provides it with the high temperature (T > 10 7 - 10 8 K) 
and the high thermal pressure that allows it to drive a strong 
shock (the so called outer shock) into the surrounding ISM 
The outer shock is responsible for sweeping and accelerat- 
ing the surrounding gas into a large-scale shell, while the 
high sound speed, thermalized ejecta, fills most of the vol- 
ume formerly occupied by the swept up gas. The two gases, 
the ejecta and the swept up ISM, are separated by a con- 
tact discontinuity and thus there is little enrichment, at this 
time, of the galaxy ISM (Tenorio-Tagle 1996). 

Several calculations (see Tenorio-Tagle & Bodenheimer 
1988, Silich & Tenorio-Tagle 1998 and references therein) 
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have shown how the growing structures may acquire a va- 
riety of shapes, depending on how the ISM is distributed. 
Fairly round and 8-shaped remnants are expected for a con- 
stant density medium and plane stratified atmospheres, re- 
spectively. Superbubbles are also known to blow out upon 
the sudden acceleration experienced when crossing super- 
sonically a large negative density gradient. At that moment, 
the accelerated section of the shell of swept up matter be- 
comes Rayleigh- Taylor unstable and fragments, while the 
hot supperbubble interior is vented into the low density sur- 
rounding gas. The escaping material, flying with its sound 
speed, would then push again the outer shock to build a new 
and even larger shell evolving into the halos of galaxies. If the 
shock reaches the outskirts of galaxies with a speed larger 
than the escape velocity of the system, it may also establish 
a galactic wind. This is in fact detected in the case of some 
nuclear starbursts (e.g. Heckman et al. 1996, Tenorio-Tagle 
& Munoz Tunon 1998). The latter event is expected to have 
drastic consequences as all the newly processed starburst el- 
ements, originally found in the hot superbubble interior, will 
be channeled out of their parent galaxy into the intergalactic 
medium. 

Superbubbles have been recognized for their large-scale 
expanding HI shells (Heiles 1979, Brinks & Bajaja 1986, 
Maschenko et al. 1999). Some of the youngest ones are de- 
tected at optical frequencies emanating from giant HII re- 
gions (Meaburn 1980, Heckman et al. 1990, Martin 1996, 
Oey 1996), and many have been recognized by the X-ray 
emitted from their hot interiors (Wang & Helfand 1991, 
Heckman 1995, Wang 1999). 

Several authors have pointed out a large discrepancy 
between theory and observations of superbubbles. Current 
predictions are based on the fact that, in the temperature 
range of ~ 10 6 — 10 7 K, the X-ray emissivity can be ap- 
proximated by a linear function of the gas metallicity. For 
the 0.1 - 2.4 keV energy band, this can be approximated by 
a constant value A x — 3£ x 10~ 23 erg cm 3 s" 1 , where £ is 
the metallicity in solar units. A simple analytic model may 
then be developed to estimate the X-ray luminosity from a 
spherically symmetric, energy dominated bubble bound by 
a cold radiative shell, and presenting a self-similar tempera- 
ture and density distributions (see Chu & Mac Low 1990 for 
a constant energy input rate, and Silich 1995 for a power- 
law energy deposition). For a constant ambient gas density, 
and constant energy input rate, the X-ray luminosity over 
the energy band 0.1 - 2.4 keV is given by 



L w = 10 36 £/(r)L 33/35 nr /35 4 9/35 



erg s 



(1) 



where L38 is the mechanical luminosity of the starburst in 
units of 10 3S erg s _1 , no is the ambient gas number density, 
and tr is the evolutionary time in units of 10 7 yr. I(r) is a 
dimensionless integral: 
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and t = Tcut/T c is the ratio of the X-ray cut-off temperature 
(or lowest limit considered for the gas temperature to lead 
to an important X-ray emission; T cut ~ 5 x 10 5 K) to the 
bubble central temperature. 

These simple considerations have been applied to a 
number of objects, from stellar wind bubbles to starburst 
galaxies. In many cases, the observed X-ray luminosities are 



in clear disagreement with the model predictions and ex- 
ceed, by much, the model predictions (e.g. Chu & Mac Low, 
1990 for LMC bubbles; Walter et al., 1998 for a superbubble 
in the nearby dwarf galaxy IC 2574: Martin & Kennicutt, 
1995 for diffuse X-ray emission from the central region of the 
starburst galaxy NGC 5253). There are several proposed so- 
lutions. Chu & Mac Low (1990) have proposed off-centered 
supernova explosions. Franco et al. (1993) included the in- 
teraction of a fragmented ejecta with the outer shell and 
Martin & Kennicutt (1995) concluded that cloud evapora- 
tion may be a dominant mechanism to increase the X-ray 
emission. 

Here we take the model one step further and show that 
the injection of energy and mass from an aging stellar clus- 
ter not only leads to the large-scale evolving superbubbles 
detected at a variety of frequencies, but also to a time- 
dependent enhanced metallicity of their hot interior. Section 
2 discusses the time dependent production of metals released 
within the hot superbubble interior. Section 3 presents the 
calculations of the energy deposited by massive coeval star- 
bursts and the implications regarding the cooling of the gas 
and the mixing of the enriched material. The resultant time- 
dependent X-ray luminosity, is given in section 4 and some 
conclusions are drawn in section 5. 



2 TIME-DEPENDENT METAL PRODUCTION 
IN STARBURSTS 

We assume that the total amount of gas ejected by super- 
novae (SNe) includes an important fraction of newly syn- 
thesized oxygen and iron, and that the extended stellar 
outer envelope has the same metallicity as the host galaxy 
ISM. We further assume that the matter violently ejected 
by SN as a starburst evolves is efficiently thermalized at a 
reverse shock and well mixed with the gas that evaporates 
from the cold radiative shell segments. However, a contact 
discontinuity inhibits its immediate mixing with the ISM, 
which is rapidly removed and accelerated by the outer shock 
(Tenorio-Tagle 1996). Consequently, the metals ejected by 
sequential SNe are to be found in the "hot cavity" or super- 
bubble interior, causing drastic changes to its metallicity. 

We assume a Salpeter stellar initial mass distribution 
(IMF) 



n(m) — fom 



(3) 



within a range of upper M up = IOOMq and lower Mi ow — 
IMq cut-off masses and a slope of a = 2.35. The normal- 
ization constant /o is determined by the total mass of the 
star cluster Msb, 



h 



(a - 2)M S i 



Mf 



Mi 



(4) 



Assuming for simplicity that massive stars loose all their 
mass as they explode as SNe (Pilyugin 1992, Pilyugin & 
Edmunds 1996), we find the total ejected mass as a function 
of the cluster age t, 



M ej (t) = /o / m 1 - a dm = Msb M " ' ' ' 



M 2 ~ a 



M„(t) 



Mf 



Mi 



,(5) 



where the mass of the stars exploding after an evolutionary 
time t, (M*(i)), has been found using Chiosi et al. (1978) 
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and Stothers (1972) approximations to the main sequence 
lifetime of massive stars: 



M*(t) = 



10(9 x 1O 6 /£) 2 M , for 30M Q < M* < 100A/ Q 
10(3 x 10 7 /i) 5/8 A^Q, for 7M Q < M* < 3OM . 



(6) 



Heavy element yields from massive stars have been con- 
sidered in a number of papers (e.g. Renzini et al. 1993). 
However, the exact values depend strongly on the adopted 
stellar evolution models. A recent attempt to incorporate 
heavy elements injection into hydrodynamical bubble mod- 
els, by D'Ercole & Brighenti (1999), has included only av- 
eraged values. We have decided to consider several possible 
scenarios. The oxygen yield as a function of the stellar-mass 
(M+) can be approximated by two different tracks: one is the 
"no wind" (NW) Pilyugin & Edmunds (1996) analytic ap- 
proximation to Maeder (1992) and Thielemann et al. (1993) 
models for the yield 



Y Q (M) = 0.094(M - 10.5) 1 



Mr: 



(7) 



which neglects the oxygen yield for stars with a mass smaller 
than 10.5 M . The second "with wind" (WW) approxima- 
tion, follows from stellar evolution models accounting for 
stellar winds (Maeder 1992, Woosley et al. 1993). In this case 
equation (|7|) can be used within the 10. 5M < A/* < 25A/ 
range, assuming a constant yield up to the upper mass limit 



M„ 



100M 



G- 



The long term iron contamination comes mainly from 
the SNIa, which produce 0.5-0.7 M of iron after the 56 Ni 
decay. However, we are interested here in the earlier stages of 
the starburst evolution and thus the iron yield from the SNIa 
has not been taken into account. The iron production from 
type II SN is highly uncertain (see Renzini et al. 1993) and, 
therefore, we have considered two extreme models: Arnett 
(1991) and Thieleman et al. (1992). In the Arnett (1991) 
model, the iron yield increases with the mass of the star, and 
can be approximated by a linear function for stars between 
10M© < M* < 4OM . 



Y Fe (M*) = 0.02 + 0.006(A/* - 10) M e 



(8) 



For other stellar masses we have assumed constant values of 
Y Fe . Y Fe = 0.02 M , for M* < 1OM . and Y Fe = 0.2 M 
for M* > 4OA'/ . In Thielemann et al. (1992), the iron yield 
is approximated by the exponential function 



Yf^M, 



0.423 



exp[0.31(Af* - 10.5)] 



0.045 Mq, 



(9) 



within the 13M - 25M range, and assumed to be con- 
stant outside this range and equal to Y Fe (M+) = 0.24 M 
for low mass stars, and Y Fe (M+) = 0.05 M for stars with 
an initial mass larger than 25 M . We shall refer to the 
above approaches for deriving the iron production rates as 
A (for Arnett 1991), and T (for Thielemann et al. 1992), 
respectively. 

We have also assumed that the stellar matter ejected 
before the first SN explosion (~ 3 Myr) has the same metal- 
licity as the ISM (see table 1). Note that the expressions 
used for the yields are for solar metallicity stars (which are 
the ones available) while the stars in the models have metal- 
licities of 0.1 Z . In what follows, the integrated metal con- 
tent of superbubbles is given in solar values for which Zo = 
0.0083 and Z Fe = 0.00126 (Grevesse, Noels & Sauval 1996). 



The metallicities of the hot bubble interior relative to 
solar are given by 

M & j t o/Zo + ZlSMM ev 



Zhot,o = 



Zhot,. 



M ev + M ej 

Mej, F e./Z F e + ZlSAlMev 

M ev + M ej : 



(10) 



(11) 



where M ev is the mass added to the shocked wind region 
due to the cold outer shell evaporation. 

Figure 111 shows the time dependent total amount of mat- 
ter injected into the superbubble interior (M e j ec t a ) as a con- 
sequence of type II supernova for the NW and A approxi- 
mations and the WW and T assumptions (panels a and b, 
respectively). This is to be compared with the total amount 
of oxygen (Mo) and iron (M Fe ) produced by the star cluster 
during the first 40 Myr of evolution. Note that at this time 
the amount of iron has not yet reached its final value, as an 
important contribution is expected later on from stars with 
a mass smaller than 10 M . 



3 SUPERBUBBLES EVOLUTION 

We have carried out calculations in order to find the im- 
pact that the metals ejected by a 10 6 M coeval starburst 
may have on the ISM of galaxies with low metallicity (Z 
~ 0.1 Z ). The calculations were carried out with our 
three dimensional Lagrangian code, based on the thin layer 
approximation (Bisnovatyi-Kogan & Silich 1995, Silich & 
Tenorio-Tagle 1998) . In all cases the total mass of the host 
galaxy amounts to 10 10 M , while the gas mass amounts 
only to 10 9 M . The galaxy model was approximated with 
a similar prescription to that used by Li & Ikeuchi (1992), 
Tomisaka & Bregman (1993), Silich & Tenorio-Tagle (1998) 
and Tenorio-Tagle et al. (1999; hereafter paper I). The gas 
density distribution allows for two isothermal components. 
One is related to a central dense molecular disk with a 
mass of 5 x 10 7 M and the second represents the low 
density, extended neutral halo. Both components are in a 
quasi-equilibrium state supported by rotation and random 
gas motions with a velocity dispersion of 20 and 80 kms -1 , 
respectively. The temperature of both components is ~ 1000 
K. 

The assumed initial density distribution of the model 
(Figure El) is similar to that of case Al in Paper I. Despite 
its large column density (~ 10 23 cm -2 ) only a small fraction 
of this mass arises from the extended halo, and most of the 
undisturbed column density is within the dense molecular 
core. The 10 6 M starburst was assumed to produce a con- 
stant mechanical energy at a rate of ~ 3.2 x 10 40 erg s _1 (see 
Leitherer & Heckman 1995 synthesis models for low metal- 
licity starbursts) for the first 40 Myr of its evolution, lead- 
ing to a giant evolving superbubble. Three different cases 
are presented. In all we assume that the supernova products 
(oxygen and iron) rapidly mix with all the gas within the 
superbubble interior. As discussed in section 2, two different 
metal production rates, for both oxygen and iron were used, 
and their possible impact is discussed here and in section 4. 

Case 1 follows the amount of mass thermally evapo- 
rated from the cold shell segments and considers its mixing 
with the ejecta from type II supernovae. The NW and A 
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Figure 2. The logarithmic density distribution for a 10 9 
M q ISM model. Labels correspond to the densities and the outer 
contour represents the edge of the galaxy. 



20 
Time (Myr) 

Figure 1. Time-dependent metal production from starbursts. 
The panels display the cumulative amount of matter (in units 
of the starburst mass Msb) ejected by SNe as a coeval starburst 
evolves in time (solid line). Also shown are the amounts corre- 
sponding to oxygen (dotted lines) and to iron (dashed lines), (a) 
The oxygen production values were derived under the NW as- 
sumption and those for iron following the A approximation, (b) 
The corresponding values assumed the WW and the T approx- 
imations for oxygen and iron, respectively. The right hand side 
ordinate indicates the mass of O or Fe in units of the starburst 



approximations were used for case 1A, and the WW and 
T approximations, for case IB. Case 2 neglects the process 
of thermal evaporation and thus the metals are allowed to 
mix only with the ejected stellar envelopes, the latter as- 
sumed to have the same metallicity as the ISM in the host 
galaxy. This calculation leads to the largest metal rich su- 
perbubbles. Case 3 is almost identical to case 1, with the 
exception that the assumed mechanical energy input rate is 



the one expected from a 10 7 M coeval starburst (3.2 x 
10 41 erg s _1 ) instead of the previous 10 6 M© ones. Case 3A 
follows the bubble evolution within a low metallicity ISM 
(Z=O.1Z0), whereas case 3B assumes solar metallicity for 
the host galaxy ISM. The model input parameters are sum- 
marized in Table pi 

The calculations follow the change of the metallicity 
in the superbubble interior with time and apply the cor- 
respondingly modified cooling rates. For the outer shell we 
assumed that its metallicity remains constant in time and 
equal to Zism- The results are compared to our case Al from 
paper I, which assumed that the cooling function, scaled 
to the initial metallicity of the host galaxy, remained un- 
changed with time, despite the obvious injection of new met- 
als into the superbubble interior. Furthermore, our compar- 
ison case Al also assumed, as in most calculations in the 
literature, that the same cooling function could be applied 
throughout the flow. 

Figures Ha and tab, show the evolution with time of the 
size and expansion velocity of the superbubble for cases 1A 
and 3A, respectively. The figure displays values of the fastest 
expansion velocity measured along the symmetry axis, as 
well as of the largest radius acquired by the superbubble dur- 
ing the first 40 Myr of evolution. The maximum expansion 
velocity shows an initial deceleration, followed by a strong 
acceleration that leads to several hundreds of kms -1 im- 
mediately after blow-out from the central disk. This short 
phase is abruptly interrupted once sufficient halo matter has 
been swept into the expanding shell causing, once again, a 
steady deceleration of the remnants. The maximum expan- 
sion speed is to be compared with the escape velocity from 
the host galaxy (V esc ), to discern whether the decelerating 
remnant remains bound (as in case 1) or is to reach the 
galaxy outer boundary and eject its contents into the intra- 
cluster medium (as in the more energetic case 3) . The results 
of cases 1, 2 and those of case Al of paper I are identical, as 
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Table 1. Model parameters 



Model 


Lss 


Yo 


Y Fc 


%ISM 


Shell evap 




erg s 






Z© 




1A 


3.2 xlO 40 


NW 


A 


0.1 


yes 


IB 


3.2 xlO 40 


WW 


T 


0.1 


yes 


2A 


3.2 xlO 40 


NW 


A 


0.1 


no 


2B 


3.2 xlO 40 


WW 


T 


0.1 


no 


3A 


3.2 xlO 41 


NW 


A 


0.1 


yes 


m 


3.2 xlO 41 


NW 


A 


1.0 


yes 


Al 


3.2 xlO 40 


- 


- 


0.1 


yes 
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all three assumed the same mechanical energy injection and 
the same galaxy. The generated remnants, in all three, have 
identical size, expansion speed and total amount of swept- 
up mass. Below we show that, on the other hand, the total 
luminosity and the metal content inside the superbubble, 
are different for the different models. 

The results for cases 1A and 2A are shown in Figure m 
The diagram shows the amount of matter swept by the outer 
shock (M s heii), and the mass of thermally evaporated gas 
from the outer shell that is ejected by supernovae (M^ot) (all 
of which is to be found in the hot superbubble interior) as 
a function of time. Also shown are the cumulative amounts 
of matter ejected by SNe (M e j ecta ) and the corresponding 
amounts of oxygen (Mo) and iron (Mf e ). The sharp M 3 heii 
decrease observed at about 5 Myr results from the outer 
shock merging at the mid-plane of the galaxy (see Silich & 
Tenorio-Tagle 1998). 

In case 1 (Figure Ma), the injected metals were allowed 
to mix thoroughly with the matter evaporated from the ex- 
panding shell. Note that for such a massive starburst, the 
amount of hot matter evaporated from the shell reaches a 
large proportion of the total mass of the shell (M^ot ~ 10% 
Mshe.il)- The total oxygen mass surpasses the 10 4 Mq value 
while the total mass in iron is ~ 10 3 M © . 

Case 2 (Figure Mb) assumes no mass evaporation from 
the shell, and thus apart from the small contribution due 
to the early wind phase, the amount of Mhot is very simi- 
lar to the matter ejected by supernovae (M e j ec ta)- Thus the 
amount of oxygen and iron are very close to the values ex- 
pected from the synthesis caused by the massive starburst. 

Figure H shows the evolution of metallicity of the super- 
bubble interior as a function of time for the different model 
assumptions regarding the yields, the mixing of heavy ele- 
ments and the mass evaporation. Case 1A (using the NW 
and A approximations) rapidly (between 4-5 Myr) reaches 
over solar metallicities Zo ~ 8 Z © for oxygen, and solar 
values for the iron tracer (Zf«0 (Figure Ha). Zo then slowly 
decays to about 3 Z © after 10 Myr of evolution, to solar af- 
ter 20 Myr and to under solar (~ 0.5 Z © ) at the end of the 
calculation. Zf b becomes under solar after 6 Myr and slowly 
approaches the metallicity of the host galaxy (0.1 Z © ). The 
metallicity that results from case IB, assuming the WW and 
T approximations for the production of metals (Figure Bp), 
reaches Zo ~ 2 Z © after 6 Myr of evolution, falling below 
solar metallicity just after 10 Myr to slowly approach 0.1 
Z © at the end of the calculation. Zp e on the other hand, 
hardly shows a significant variation in this case and remains 
under solar throughout the calculation. 

Case 2 (without mass evaporation) reaches, within short 




10 20 30 40 

Time (Myr) 

Figure 3. Evolution of the superbubbles. (a) Case 1A and (b) 
case 3A; shown as a function of time are the maximum expansion 
speed (kms -1 ) measured along the symmetry axis (solid lines), 
the galaxy escape velocity (dash-dotted lines) and the radius (in 
kpc) reached by the superbubble. 
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Figure 4. The log of the swept up mass (in Mq) and that inside 
the superbubbles for cases 1A (a) and 2A (b), as a function of 
time. The curves represent the log of: total swept up mass (solid 
lines), total mass inside the superbubbles (dotted lines labeled 
Mhot)- The latter includes both the amount of matter thermally 
evaporated from the shell (not shown) and the mass ejected by 
SN (dashed lines labeled M e j ec t a , coinciding with M^ot in panel 
b). Also shown are the log of the total oxygen and iron mass 
(dash-dotted lines). In both cases these were calculated following 
the NW and A approximations (see section 2) 



times after the start of the SN phase (3 Myr) , a value of Zo 
larger than 30 Zq (case 2A). Even after a substantial fall, 
as the less massive stars become SN, this case ends with 
Zo ~ 15 Zq (Figure nc). In this case Z^e ~ 6 Zq after 
6 Myr of evolution and remains almost constant throughout 
the calculation. Case 2B (under the WW and T assumptions 
for the metal production rate) presents a rapid rise in Zo 
reaching a maximum Zo =8Zq at about 10 Myr and then 
slowly declines to 6 Z q after 40 Myr. The Z^e on the other 



hand, remains subsolar for the first 12 Myr and then rapidly 
rises to reach 6Zq after 40 Myr (Figure pp). 

For the most energetic starburst (case 3) large values 
of Zo > 10 Z © are found within the first 15 Myr (Fig- 
ure He), when the NW and A metal production rates and 
Z;sm = OTZq are used (model 3A). The enhanced metal- 
licities however, are not diluted in this case as the fast and 
hot quasi-adiabatic outer shell has inhibited mass evapora- 
tion. 

As we increase the initial ISM metallicity from 0.1 Zq to 
Zq for model 3B, the interior bubble metallicity is lower and 
more peaked with time (Figure B). This unexpected result 
is produced by a much faster shell cooling, as it is more 
metal rich, and thus the sarburst ejected metals are more 
efficiently diluted by the evaporated gas from the shell. 

Finally, we comment on present day abundance deter- 
mination of X-ray gas. ROSAT observations allow, in prin- 
ciple, for the derivation of metal abundance, of a (T ~ 
1 — 30 x 10 6 K) plasma, by X-ray spectral fittings. How- 
ever, the process is problematic, because of the poor en- 
ergy resolution of the ROSAT-PSPC detector, and the de- 
rived ISM metallicities disagree, in many cases, with those 
deduced from optical observations (Trinchieri et al. 1994; 
Saracco & Ciliegi 1995; Bauer & Bregman 1996). The sit- 
uation has somewhat improved with the ASCA and Bep- 
poSAX observations yet the X-ray measurements still sug- 
gest significantly sub-solar abundance for Fe combined with 
somewhat sub-solar abundances for a-elements (Bauer & 
Bregman 1996, Ptak et al. 1999, Persic et al. 1998). In some 
of the cases, the optical observations clearly indicate higher 
metallicities. 

In other words, the complexity of the X-ray spectra of 
starburst galaxies, combined with the poor spectral resolu- 
tion of the pre-Chandra and pre-XMM observations, makes 
the determination of metallicities very model dependent 
(e.g. Netzer 1999). Dahlem, Weaver & Heckman (1998) and 
Weaver et al. (1999) have discussed the difficulties in deriv- 
ing metallicities from X-ray data, and, by combining ROSAT 
and ASCA data for nearby edge-on starbursts like NGC 253 
and M 82, they concluded that the fitting is consistent with 
near solar abundances and therefore, extremely low metal- 
licities as derived by only ROSAT or ASCA data are no 
longer required. 



4 THE X-RAY LUMINOSITY OF 
SUPERBUBBLES 

The X-ray emission from superbubbles arise from two differ- 
ent physical regions: the dense outer shell of accelerated ISM 
and the hot superbubble interior (Suchkov et al. 1994, Silich 
& Tenorio-Tagle 1998, Strickland & Stevens 1998, D'Ercole 
& Brighenti 1999). Our description of the inner bubble struc- 
ture is based on the Weaver et al. (1977) similarity solution, 
and assumes that the density n and temperature T profiles 
can be approximated by 



n = n c (l - 

t = r c (i 



(12) 
(13) 



where n c and T c are the central hot gas number density and 
temperature, and x = r/R 3 h is the dimensionless distance 



© 0000 RAS, MNRAS 000, 000-000 



X-rays from superbubbles 7 




Li 



10 20 30 

Time (Myr) 




Figure 5. The panels display the oxygen and iron content of superbubbles (in solar units: Zo = 0.0083 and Zp e =0.00126, from Grevesse 
and Noels 1996) for the models (a) : 1A (with thermal evaporation of the cold outer shell segments) and NW and A approximations; (b) 
IB with WW and T approximations; (c) 2A (no mass evaporation, NW and A approximations); (d) 2B, (no mass evaporation, WW 
and T approximations); (e) 3A, a 10 7 Mq starburst within a O.IZq ISM, and (f) 3B, a 1O 7 M0 starburst within a Zq ISM. 
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from the bubble center. To account for the expected differ- 
ence in density and temperature distributions for adiabatic 
and radiative shell segments, we approximate the inner bub- 
ble luminosity as 



eL R + (1 - e)L, 



(14) 



where e is the ratio of the radiative surface segments to the 
total remnant surface area, and L_r and La are the X-ray 
emissions from the interior of the fully radiative and adia- 
batic bubbles, respectively (Silich & Tenorio-Tagle, 1998) 



L xM = J^ &<t>$* smOdO J Rcut £n 2 (r)A x (T)r 2 dr 
3£A-Vr2/(T C ,A) 

Here I(T C , A) is a dimensionless integral given by 



\ T c J 



dT, 



(15) 



(16) 



Kcut and Tcut are the X-ray cut-off radius and temperature 
respectively. We used self-similar power index A = 2/5 for 
Lr, and A = 1/20 for a more homogeneous gas distribu- 
tion within a bubble with an adiabatic hot shell. The X-ray 
luminosity from the adiabatic shell segments is taken to be 



shell — £ y^n shock A x (T 3 t loc k)ARd'£, 



(17) 



where dS is the adiabatic segment surface area, and AR is 
the segment thickness. 

Figure Br shows the total and the shell contribution to 
the X-ray emissivity of case 1A, using the hot gas metallici- 
ties shown in Figure Ha. The contribution of the shell to the 
total emission becomes most important immediately after 
blow-out of the remnant from the central gas distribution 
into the extended halo of the host galaxy. During blow-out 
the shock reaches speeds of several hundreds of kms -1 lead- 
ing to the high temperatures (T ~ 1.4 x 10 7 v 2 hock ; where 
Vshock is the shock velocity in units of 10 3 kms -1 ) which al- 
low the shocked gas to radiate in the X-ray regime. The shell 
contribution becomes rapidly less important as the shock 
slows down and the shell cools below the X-ray cut-off tem- 
perature (~5x 10 5 ) and condenses into a narrow and dense 
outer boundary of the superbubble. On the other hand, the 
luminosity of the superbubble interior is found to remain 
within the range of observed values (10 38 - 5 x 10 39 erg 
s _1 ) throughout the evolution. Figure nb compares the time 
dependent X-ray luminosity for our case 1 under the two 
assumptions for metal production rate. Clearly, the lower 
metallicities reached in the WW model, lead to a smaller X- 
ray emission, although the difference never amounts to more 
than a factor of 5. 

Figure He gives the contribution to the X-ray luminosity 
of case 1A in two different energy bands, as a function of 
time. The two energy bands considered are from 0.1 to 2.2, 
and from 1.6 to 8.3 keV (X-ray emissivities from Suchkov et 
al. 1994). Clearly, most of the X-ray luminosity arises from 
the soft X-ray component. The high energy X-rays are due 
mostly to the larger central temperatures in the superbubble 
interior. 

Figure m shows the time evolution of Lx , per unit mass, 
in models 1A, 2A, 3A and 3B, as well as the same quantity 
for model Al from paper I which did not include changes in 




10 20 30 40 

Time (Myr) 

Figure 6. The log X-ray luminosity, as a function of time, for 
case 1A under NW and A approximations, (a) shows the outer 
shell contribution (dotted line) and the total luminosity (solid 
line); (b) compares the above total luminosity with that derived 
for case IB under the WW and T assumptions (dotted line); and 
(c) shows the X-ray luminosity in the two indicated X-ray bands 
for case 1A. 



© 0000 RAS, MNRAS 000, 000-000 



X-rays from superbubbles 9 



' 32- 



..v"v.-1 






J 


34 




.' ' V: ^^ 




1A 


' i: 






t 
t 




Al 


s 
\ 
\ 

" - 2A 

i . . . . i . . ~T -. , i 



10 20 30 40 

Time (Myr) 

Figure 7. The log of the X-ray luminosity of superbubbles per 
unit mass as a function of time for cases 1A (low Z), comparison 
Al (from paper I) (low Z), 2A (low Z), and massive starburst 3A 
(low Z) and 3B (solar abundance), as labelled. 



metallicity due to the injection of new metals into the cavity. 
The main difference between case 1A and reference case Al 
comes from the hot supcrbubble interior and can amount to 
more than an order of magnitude, particularly after 10 Myr, 
once the shell contribution has become negligible in both 
cases. Note that the X-ray luminosity of case Al is in fact 
rather similar to that of case 2A (without thermal evapora- 
tion). However, in case 2 the reduced X-ray emissivity from 
the interior is due to the fact that, without the substantial 
extra mass input, the amount of matter able to radiate is 
orders of magnitude too small. Case 3, our most powerful 
starburst, leads to the largest luminosities of all cases. 

It is interesting to notice an important difference be- 
tween the bubble evolution in low (case 3A) and high (case 
3B) metallicity ISMs. Up to 10 Myr, the X-ray luminosity 
of the bubble in the metal-rich galaxy is almost an order 
of magnitude larger than that for the low metallicity one. 
However, as the high metallicity bubble reaches the radia- 
tive phase much earlier than the low metallicity one, the 
metal-rich shell (3B) cools rapidly around 11 Myr, whereas 
in the metal-poor case (3A) it remains adiabatic and hot. 
This results in a rapid decrease of ~L X in the high metallicity 
model, which, after 12 Myr, drops below the X-ray emission 
from the low metallicity one, and decline more slowly af- 
terwards. The plot indicates that similar X-ray luminosities 
per unit mass (~ 10 33 erg s _1 Mq ~ 1 ) arise from cases that 
allow mass evaporation from the outer shell, irrespective of 
the mass of the starburst causing the superbubble. This is 
very different from the calculations that do not account for 
the change in metallicity of the superbubble interior (case 
Al, from paper I) and also from those that do not account 
for mass evaporation (case 2). 



5 CONCLUSIONS 

The calculations presented here, with the exception of those 
for the more massive starburst (case 3), differ only slightly 
in their final outcome regarding the size, expansion veloc- 
ity and amount of matter swept by evolving superbubbles. 
However, they show enormous differences in the metallicity 
of the hot interior, and thus produce largely different X-ray 
emissivities. 

The effects of injection of new metals into the superbub- 
ble interior is most noticeable during the HII region lifetime 
(the first 10 Myr). If oxygen is used as a tracer, the metal- 
licity reaches, immediately after the start of the SN phase, 
values well above solar metallicity. The maximum value is 
achieved 3 to 5 Myr after the beginning of this phase. Later, 
mixing with the evaporated mass slowly reduces the impact 
of SNe. Note however, that if massive stars are assumed not 
to have strong winds, large values of Z > ZqwouFI remain 
present inside the superbubble for at least the time required 
to reach the end of the type II SN phase (40 Myr). This is 
true even in the extreme case where the host galaxy initial 
metal abundance is well below solar. On the other hand, 
if iron is used as the tracer of metallicity, values of Z^e < 
Z q are always predicted, regardless of the assumed iron pro- 
duction scheme. The above results are to be compared to the 
cases without mass evaporation, all of which led to metallic- 
ities larger than solar throughout the evolution, regardless 
of the assumed tracer. 

More massive starbursts result in very strong shocks 
and longlasting hot, quasi-adiabatic outer shell from which 
thermal evaporation is strongly inhibited. Such starbursts 
show over solar metallicity throughout their evolution, re- 
gardless of the assumed tracer. 

The most important conclusion is that enhanced metal- 
licity of superbubble interiors can strongly influence their 
X-ray luminosity, bringing theory and observations into a 
better agreement. 
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